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We rederive the thermodynamical properties of a non interacting gas in the presence of a minimal 
uncertainty in length. Apart from the phase space measure which is modified due to a change of 
the Heisenberg uncertainty relations, the presence of an ultraviolet cut-off plays a tremendous role. 

The theory admits an intrinsic temperature above which the fermion contribution to energy 
density, pressure and entropy is negligible. 

I. INTRODUCTION 

Modern physics is an edifice in which every stone is tightly linked to the others. A slight modification in one area 
may produce important changes in different fields. Quantum mechanics starts with the commutation relations. Once 
they have been fixed, using a representation of the operator algebra we can in principle solve the Schrodinger equation 
whose solutions potentially contain all the physics of the (non relativistic) system studied. The correspondence 
principle which states the link between the commutator of two variables of the phase space and their classical Poisson 
bracket is one of the basic axioms of quantum mechanics. It is not deduced from another assumption and cannot be 
judged in an isolated manner: only the whole theory can be confronted to experience through its predictions. 

Some authors have investigated the consequences of the alterations of these commutation relations on the observables 
of some physical systems. In particular, some deformations of the canonical commutators studied by Kempf-Mangano- 
Mann(K.M.M) induce a minimal uncertainty in position(or momentum) in a very simple way, providing a toy model 
with manifest non locality [1]. The implications of some of these quantum structures have been studied for the 
harmonic oscillator [1-3] and the hydrogen atom [4]. The transplanckian problem occurring in the usual description 
of the Hawking mechanism of black hole evaporation has also been addressed in this framework, for the Schwarzschild 
and the Banados-Teitelboim-Zanelli(BTZ) solutions [5,6]. These studies have extended the work on the entropy of 
the black hole [7] and Hawking radiation [8] in theories with modified dispersion relations [9-13] 

One of the purposes of this article is a re-investigation of the modifications these models induce, not in the charac- 
teristics of a single particle, but in the behavior of a macroscopic system. The modern presentation of thermodynamics 
basically relies on statistical physics. According to the system under study(isolated, closed or opened), one uses an 
ensemble (micro canonical, canonical or grand canonical) and the corresponding potential(entropy, free energy, grand 
potential) to derive thermodynamic quantities(pressure, specific heat, chemical potential, etc). 

The thermodynamic potentials are related to the derivatives of the partition function which itself is an average(of 
a quantity which depends on the ensemble used) on the phase space. To define a measure on the phase space, one 
needs to know the extension of the fundamentals cells. In the "classical" statistical physics, the Heisenberg uncertainty 
relations are used to show that this volume is the cube of the Planck constant and the undiscernability of particles 
justifies the Gibbs factor. This can be inferred from quantum statistical physics in which the sum giving the partition 
function is made on a discrete set of states. 

If one modifies the commutators, one changes the Heisenberg uncertainty relations. The measure on the phase space 
is no more the same; this results in new partition functions and consequently different thermodynamical behaviors. 
From the quantum point of view, the energy spectrum of systems are modified by the change in the commutation 
relations. 

The thermodynamics of some models displaying a minimal uncertainty in length has been analyzed before [14-16]. 
We shall give here better approximations and correct some sign errors. The implications for the early universe have 
been at the center of many works among which [17,18]. The standard big bang in a universe in which an ultraviolet cut 
off appears in a toy model exhibiting a modified dispersion relation has similarly been analyzed [19]. The difference 
between the two approaches relies in the fact that the dispersion relations are different; in the first models, they come 
from an assumption on the structure of the quantum phase space. 

In these works, the equation of state used for radiation was obtained considering bosons. In the usual case the 
contributions of fermions to energy, pressure and entropy are simply the seven eights of the ones corresponding to 
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bosons. This renders the equation of state insensitive to the ratio between fermionic and bosonic degrees of freedom. 
We show this to be drastically changed in the new framework. 

Many studies have been devoted to cosmological perturbations in transplanckian physics [20-31]. The consid- 
erations we develop here may, in this context, be seen as relevant only in the pre-inflationary era. However, the 
phenomenological bounds [4] are much lower than the Planck scale. If one adopts a less restrictive point of view, then 
the new scale may generate some seizable effects. 

The article is organized as follows. In the second section we give a very brief survey of the two models we will 
study. They possess a minimal length uncertainty and so the quasi-position representation plays a crucial role. The 
third section is devoted to the study of the non interacting gas. We obtain its equation of state, entropy, etc. We use 
its specific heat to set a bound on the deformation parameter and find it to be in agreement with previous ones. In 
the fourth section we address the question of the black-body radiation in the new context and find sensible differences 
between fermions and bosons at very high temperatures, like for the ultra relativistic gas. 



II. DEFORMED THEORIES 



As shown by KMM [1], the quantum mechanical theory defined by the commutation relation 

[±,p] =ih(i + pp 2 ) (l) 

is endowed with a minimal length uncertainty Ax m i n — h\f]5. The presence of this minimal length uncertainty 
implies that no position representation exists. The concept which proves to be the closest to it is the quasi-position 
representation in which the operators are non local: 

p = tan (ihy/^d^j , x = £ + h^/~3 tan (ihyffid^ . (2) 

The most obvious extension to a three dimensional space is obtained by taking the tensorial product of three 
such copies. It will be referred to as the A\ model. It has translational invariance but lacks the symmetry under 
rotations. Another extension preserves rotational and translational symmetry; it will be referred to as the A 2 model. 
Its commutation relations are 

[x j ,p k ]=ih{f(p 2 )6 j , k +g(f)p jPk ) , g(p 2 ) = P , f( P 2 ) = - ^ = • (3) 

We will need the form of the momentum operators in this model: 

M = -,ng(— A] - , where A = £^ . (4) 



III. THE NON INTERACTING GAS IN DEFORMED THEORIES 



Introducing the momentum scale f3 , it is straightforward that with the Boltzmann constant k , the light velocity c 
and the mass m of any particle, one can construct on purely dimensional grounds the characteristic temperatures 

1 c 

The first is "non relativistic" , particle dependent while the second is "relativistic" (c dependent) and universal. We 
are interested in what happens at and above these temperatures. 

The constant (3 which appear in the KMM algebra is a free parameter. What numerical value can it assume? It 
was suggested [32] that the minimal length uncertainty h^f]3 should be of the order of the Planck scale. We shall 
adopt a less restrictive point of view here. The only constraint is that the deformed commutator should not lead to 
contradictions with the predictions of the orthodox theory which have been observed experimentally Our attitude 
is inspired by recent works which have shown that a new physics may take place well before the Planck scale [33]. 
Taking the most stringent phenomenological constraint hy/f3 < 10~ 16 m [4], one finds a lower bound to the characteristi 
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crelativistic temperature : T cr > 10 13 K . If one assumes the minimal length uncertainty to be of the order of the 
Planck scale, then T cr = 10 16 K. 

We now briefly summarize the formula of statistical physics we will need. In the usual theory, a system which 
is in contact with a large heat reservoir and doesn't exchange particles with the surroundings has to be studied in 
the canonical ensemble [34], [35]. Its equilibrium state will be described by a fixed temperature and a fixed particle 
number while its energy will fluctuate around a mean value. Strictly speaking, for such a system the particle number 
N is fixed once and for all. But, one knows that when phase transitions are not present, the description given by the 
canonical and the grand canonical ensembles are very close. This will be used to compute the chemical potential in 
the canonical ensemble where the calculations are easier. 

The most important quantity will be the canonical partition function Z(T, V, N) which is defined in terms of the 
Hamiltonian operator H by the equation 

Z(T,V,N) = Trexp(^-^j . (6) 

The free energy is related to the partition function by 

F(T,V,N) = -kTlnZ . (7) 

In these variables the pressure P, the entropy S, the chemical potential [i, the internal energy U and the constant 
volume specific heat Cy read 

ap fjp ftp fjrr 

P = ~ W , S = -- , , U = F + TS , and Cy = - . (8) 



A. The A\ model 

1. Low temperatures 

Low temperatures correspond to non relativistic behaviours. Let us study a non interacting non relativistic gas. 
One has to solve the Schrodinger equation for a particle in a cubic box of length L. This will be done in the quasi 
position representation because of the lack of a position one. For simplicity, let us first consider a one dimensional 
system; one obtains the solution 

i>(t, ~ e' lEt exp (±^j arctan V 2m P hE ^j ■ (9) 
When the periodic boundary condition i[)(t,€ = 0) = ^>(t, £ = L) is imposed, one finds the quantification of energy 

E ^2r^ t ^{-^—) ' (10) 

(n being an integer) already obtained in [1,5]. This leads to a cut off in order to avoid a non monotonic dispersion 
relation: 



l>sup F 



(11) 



E being the integer part function(not to be confused with the energy). 

This cut off n sup is necessary in order to prevent a divergence of the partition function which would take place 
otherwise, due to the periodic nature of the energy(Eq.(10)). But this still allows an infinite energy provided that 
the length of the box is fine tuned in such a way that the result of its division by the minimal length uncertainty is 
an integer. This would not be problematic since such an energy would have a vanishing contribution to the partition 
function. This conclusion will remain untouched at high temperatures. 

The one particle partition function is given by the formula 



n=0 n=0 



(12) 
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The sum on n can be approximated by an integral on dn if the size of the box L is big enough. Introducing the 
integration variable p by 

n= ^^ arCtan (^) (13) 
and replacing, as usual in statistical physics, the length L by an integral on position, we find 

Z{T,V,l) = \J dxdp Y~tfitf- e ~ p2/2mkT ■ ^ 

The formula given in Eq.(14) admits a semi classical interpretation. Let us first consider its limiting case (3 = 0. 
Classically, the system can be seen as a point evolving in phase space. The probability for the system to be in 
a configuration in which the first particle is in the region q\ ± dq\,p[ ± dpi , the second particle in the region 
02 ± dq 2 ,p2 ± dp 2 , • • • is proportional to e~ E ^'^ and proportional to the number of elementary cells contained in 
the volume of the aforementioned region. At the quantum level, the Heisenberg uncertainty relation of the usual 
thcory(written in the one dimensional case) ApiAqi > Ti/2 assigns to each elementary cell a volume h. The number 
of such cells contained in the region under consideration is therefore 

n V d 3 p t d 3 qi 
—f^~ ■ (15) 

2 — 1 

One sees that in the new theory, one can simply keep the usual dispersion relation and modify the elementary cell 
volume. At the quantum level, this appears as a Jacobian linked to the change of variables(n — > p). This could be 
anticipated with a semi classical reasoning. The new Heisenberg uncertainty relation implies 

AxAp> |(l + /3(p 2 )) . (16) 

It assigns to the elementary cells of the phase space of the new theory a volume h(l + j3p 2 ) which replaces the usual 
factor h. From this we could conjecture a simple recipe when dealing with the semi classical approximation: it is 
obtained by keeping the classical dispersion relation but modifying the measure in a way consistent with the new 
Heisenberg uncertainty relation. 

However, it should be noted that, in the new theory, the range over which one integrates is finite, due to the presence 
of the cut-off which depends on the volume but not the temperature: 

Due to the equality 

lim« = l , (18) 

the upper bound goes to infinity and the volume of the elementary cell tends to the usual one as the deformation 
parameter is sent to zero. It can be anticipated that, because the integrand of Eq.(14) is a rapidly decreasing function, 
taking the upper bound to be infinite will not introduce an appreciable error in most cases. 

Going to the three dimensional extension A\ , one has simply to take the product of the elementary cells in the 
three directions. The integral over the positions is obvious; the change of variable w = p 2 /2mkT, allows one to write 
the one point partition function as 



Z(TA'. \ ) = ^J , with ./ i / l ' x l> : ~ : ~ 1 2 ( 1 ' -y--) d(jJ \ 
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A= GirJ ' ^ = YY ■ (19) 

When (3 vanishes, the upper bound equals infinity and the integral giving J assumes the value one so that the 
undcformed theory is recovered. The total partition function Z is found to be related to the usual one (corresponding 
to /? = and now denoted Z*) by the relation 
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The free energy then becomes 

F = — NkT In J . 

The thermodynamical quantities are affected in the following way : 



(20) 
(21) 



P = P* + NkT 



U = U*+ NkT' 



l_dJ 
1W 
ldJ_ 
Jdf 



S = S* + NklnJ + NkT 



id£ 

Jdf 



H = m* - fcTln J - NkT 



1 dJ 
J ON 



Cy = Cy 



1 dl 

2NkT— — + NkT 2 



1 



dJ 



J 2 \dT 



1 d 2 J 
JdT 2 



(22) 



As can be seen from Eq.(19), J depends on the volume only trough the cut-off whose influence will be seen to be 
negligible for the system under study. Thus, the equation of state PV = NkT will remain valid, thanks to Eqs.(22). 
The presence of the temperature and the absence of the number of particles in the expression of J results in the fact 
that the entropy receives two contributions while the last term of the chemical potential in Eqs.(22) vanishes. The 
internal energy is also modified and by way of consequence the specific heat at constant volume too. 

The integral giving J can not be computed analytically. However, the qualitative features of the theory can be 
obtained quite easily. When the deformation parameter goes to zero, J assumes the value one. As we shall show in a 
moment, the following parameterization holds: 



+ 



Working to second order, we find 

TdJ T , 
Jdf= ai T r + ( 2a2 



log 



ai- 



2a 2 



The last formula's interest lies in the fact that it allows a more compact expression of the entropy: 



S = Nk^- + \oga j +Nk\og j- 
so that an adiabatic process takes the form: 



V {2irmkT\ 



3/2 



{—) 



T 

1 + 2<7l — 



3cr 2 + -0i 



The formulas displayed in Eqs.(22) are used to recast the equation of state and the specific heat as 



2 T, 
Considering the reaction 



3 T , 2N (T 



Cv = Nk 



^+201- + (6a 2 -3a 2 ) — 



T 



aA + bB^dD + eE 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



the expression of the chemical potential /z shows that the densities at equilibrium Xi = Ni/V obey the modified law 
of action of masses 



yd ve 
A D A E 

A A A B 



c tej.3/2(-o-6+d+e) 



T T T 

1 + 2di7i ) £) — h 2e<j\,E — 2aai t A 7^ — 



lc,D 



Lc.E 



Lc,A 



2ba X) i 



lc,B 



(29) 



since each particle, having its own mass, possesses a specific critical temperature. 

Let us now evaluate the integral J in order to have numerical estimates of the constants a. We shall use the Mac 
Laurin development with a remainder. Let us separate the integral of Eq.(19)in two parts: 
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J3 = h + I 2 = — 



+ 



On the second interval, the inequality 



1 1 



can be used to obtain the majorization 
The asymptotic formula 

r(a,z) ~ z^e^ for 

shows that 



1 Tc_ 

2' 2T 



oo 



2T 



-1/2 



exp 



2T 



■ c st u su p exp (-uJsup) ■ 



(30) 



(31) 



(32) 



(33) 



(34) 



Since, for a reasonable /?, the characteristic temperature T c is reasonably expected to be very high, the integral I2 
which contains the influence of the upper bound can therefore be neglected. 

To evaluate I\ , we shall use the Mac Laurin theorem which states that for any sufficiently regular function / defined 
on an interval [0, a] and for any point uo on that interval, there exists another point 6(u) on the same interval such 
that 



f(u) = /(o) + /'(o)w + if (o)c 2 + ir(tfHV 3 



This gives 



This enables us to find 



{T, 
l: 



2T 



l + 2jrUJ 



2T 



1 



duie w lo 



«,.,-l/2 



2T 



2T 



< max 



On the interval involved, the following inequality is verified: 



I max f"(x) I < c s 



2T 



(35) 



(36) 



(37) 



(38) 



The integrals remaining in Eq.(37) can be expressed in terms of complete and incomplete T functions. The final result 
reads: 



1 



3\ 2T 



h = -= r - -r - — + t - — 



2 T c . 



2T\ 



(39) 



From this one reads the values of the coefficients a . 

The predictions of usual statistical physics are known to be quite accurate in ordinary conditions. We shall use this 
to give an estimate of the bound thermodynamics imposes on the deformation parameters. Consider the Helium whose 
specific heat at constant volume assumes an experimental value comprised between 12.39 and 12AUK ~ 1 mole~ 1 . The 
undeformed theory assigns the value 12. 47 JK~ 1 mole~ 1 to any non relativistic gas. Thanks to Eqs.(27), we can write 
the specific heat in the new theory as 
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C v = 12.47 (l + <7ij) ■ (40) 

The measured value tells us that 12.39 < Cv < 12.41. Assuming T = 300K , the prediction of the model does not get 
out of the experimental bounds provided that T c > 10 6 K which induces (3 < 10~ 45 . One then finds a minimal length 
uncertainty 7 < 10~ 12 meters which does not disagree with the bound derived from atomic physics considerations 
7 < 10~ 16 m [4] but is less precise. It should be stressed that this only gives an idea of the order of magnitude since 
we did not include interactions between the atoms. The estimated temperature T c at which something new should 
happen is too high for the non rclativistic approach to be reliable. Thus the interest of the next subsection. 



2. High Temperatures 



The deformed Klein-Gordon wave equation for a massless particle in a box leads to the spectrum 



p 2 — 



tan 2 ( TiV/3 



2irn x 



tan 2 (fi^ffi 



(41) 



Replacing the sum by an integral, one obtains for the factor J appearing in the partition function the expression 

r^sup r^sup r^sup 

J = dx dy dz ip(x, y, z) , (42) 

Jo Jo Jo 



where 



Viz, y,z) = i { (l + ^f) (l + ^) (l + ^f) (cxp ( ^TWT 



, uj sup = —t . (43) 



Let us first consider a temperature verifying T/T cr < 1, with an upper bound which is practically infinite 
(T cr t/T > 1, t being given by Eq.(17)). In this case, replacing the domain of integration, a cube, by a sphere should 
not introduce an important error. The result reads 



J = 1 - 12 



C(5) / T 



C(3) \T t 



C(7) / T 



7^- + — 



C(3) \T C 



(44) 



From Eqs.(22), one finds the equation of state and the expression of the entropy 



C(3) \T C 



288 



4 c(7) _ (cm) 
. c(3) U(3)y 



T 



S = ANk + NkXoglSir- 



V fkTY 



N \hc J 



1 - 36 



C(5) / T 



C(3) \T t 



288 5 



.07) /C(5) 
'C(3) + U(3) 



T 

c r 



(45) 



Like in Eqs(25,27), one obtains small departures from the unmodified theory. 



3. Very high temperatures 

What happens at very high temperatures? Like in the preceding subsection, the most salient features can be 
captured from the behavior of J. As the temperature is increased, the form of the its integrand and its upper bound 
(Eqs.(43)) show that J goes to zero. An approximation of the form 



J = CT„ 



T 

J cr 



71+1 



+ C„ +2 



T 



n+2 



+ 



(46) 



with n a positive integer will hold. Keeping the first two corrections one shows, by computations similar to those of 
the preceding subsection, that the entropy takes the form 



S = ANk + Nk log <^ 8ire- n cr, 



k 3 V 

rpn rp3 — n 

(he) 3 cr N 



1+ 



1 a. 



n+l 



>n+2 



T 



(47) 
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(from which the equation of an adiabatic process can be deduced) while the equation of state reads 



(3 - n) - — — + 



The law of action of masses now assumes the form: 



7 «+i _ o °"»+ 2 



L cr 



XpXg _ ^st rp(3-n)(d+e-a-b) 



l + (d + e-a-&)(±^±i-^ 
2 tr cr n 



J cr 



(48) 



(49) 



Let us find the integer n. One has T/T cr < 1 and T cr t/T < 1. Now, the domain of integration is small and so 
the exponential appearing in the function ip can be expanded in polynomials. The function 1 + T 2 r /T 2 x 2 admits two 
Taylor expansions; the region in which x < T/T cr will be denoted A; the other one will be denoted B. The same 
situation occurs for y and z; this leads to a partition of the domain of integration. The most important contribution 
reads 



J AAA 



cr 
T 3 



(50) 



The dependence on the volume is negligible. Comparing with Eq.(46), one has n = 3 so that the equation of state is, 
to first order, p ~ 0. 

In this theory, statistics will play a role at high temperature. As is evident from Eq. (42,43), J will go to zero as 
the temperature increases. The Bose-Einstein distribution 



9i 



1 



exp - v) 

reduces to the Maxwell-Boltzmann's one only in the limiting case 

Thanks to the relation giving the total number of particles N = ^ N i7 one is led to the condition 



N 



N 



e = 



Z{T,V,l) 



< 1 — < 8?r — 



V 



kT 



he 



(51) 



(52) 



(53) 



One concludes that neglecting statistics is accurate, at very high temperatures, only for systems whose densities are 
very small (J ~ T~ 3 ). If this is not the case, the appropriate integrand in the evaluation of J for bosons for example 
is 



(54) 

The normalization constant m ensures that J = 1 in the undeformed theory. One finds the dominant part is given by 

2 



J aaa = (- \ - log(2 V2) + 3 log(l + V3)) 



(55) 



In the formula giving J, one now has n — 2 and the associate equation of state for a gas of bosons takes the form 
p = P. 

We treat in more detail the A2 model in the following section. Our choice is mostly due to the fact that the A2 
model, possessing spherical symmetry, gives J as a one dimensional integral, contrary to the A x model(see Eq.(54)). 
Apart from leading to simpler formulas, this model seems also more suitable to the treatment of a Robertson- Walker 
universe because of this rotational symmetry. 
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B. the A2 model 



In this case, the action of the first position operator on the plane wave tpk(t, £) = exp (—iEt + fci£i + k 2 ^2 + ^£3) 
reduces to 



_ 00 /fi 2 ^ \ 

P*(t,0 = -*(«)E y* 2 

r=0 ^ ' 



(56) 



The sum of the left hand side converges only when h 2 k 2 < 2/\fj3; this is the cut-off. As shown in the last subsection, 
we do not learn much from the deformed non relativistic theory; we then go directly to high temperatures. The 
solution to the wave equation gives the dispersion relation 



E = chl ( 1 + 



(57) 



from which one infers the quantity controlling the departure from the unmodified theory, for fermions and bosons: 



J bo 



1 



2C(3) Jo 



72% 



dxx 



exp 



1 + \ 



Jfe 



3C(3) Jo 



dxx 



exp 



1 + l^-x 2 



(58) 



For temperatures smaller than T cr , a Taylor expansion of the term under parentheses and an approximation of the 
upper bound of the integral by infinity holds. This leads to the following expressions for the equation of state and 
the entropy in the bosonic case: 



C(5) / T 
= 3 + 60-^-( — 



bo 



C(3) \T 



360 21 



V fkTY 



C(7) 
C(3) 



05) 
C(3) 



S b o = 4Nk + Nklogi8n- i^—j 



C(5) 



1 + 90^ — 



C(3) \T 



T 

- 1 - rr 



+ 450 f 21^ + 4 



C(3) 



C(5)V 
C(3)J 



T 

-L r.r 



(59) 



For fermions, the behavior is roughly similar but the details are different: 

H 441-%+ 

2 ^ C(3) 

3 

N V Tic 



_ U7 ,((5) / T 



S ie =^k + Nk\oA%v- — 



1 



- 125 
225 C(5 



/C(5) 
U(3) 



2 C(3) V^c 



T 

-L r.r 



225 / C(7) 
— 441> V ' 



C(3) 



100 



C(5) 
C(3) 



T 



(60) 



When the temperature becomes of the order of the critical temperature, an important change takes place. The 
domain of integration for J in Eq.(58) becomes small and so one can approximate the integrand by its Taylor expansion 
near the origin. The difference between bosons and fermions enters into play through the difference of signs ± which 
leads to different powers in terms of the temperature. Developing the full integrand in Eq.(58) to the fourth order in 
x, one finds 



Jbo = 



1 



2C(3) 



Tcr_ 

T 



T 



Jfe 



3C(3) 



1 (Tcr 

6 V T 



I T 



16 V T 



(61) 



A computation to an order greater than the one to which we have limited ourselves brings in small corrections to the 
coefficients 13/16 , etc. Using Eqs.(48), one then finds the expression of the entropy and the equation of state. The 
difference is significant between the two statistics as can be seen from the following equations: 



P 

VJbv 



2 /—T rr 8 / T C r 

1 + -V2— + — — 

9 T 81 ^ 



P 

VJ fe 



3 T cr 
8~ 



9 

64 



Tcr 

T 



(62) 



A numerical computation supports our approximation scheme. 
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12 3 4 

FIG. 1. - plotted as a function of for bosons and fermions 

In Fig.l we plot the ratios between the density and the pressure for a gas of bosons (fermions). For temperatures 
smaller than the critical one, these ratios are close to their known value 3 in the undeformed theory. The two functions 
then rise above this value as predicted by Eqs. (59,60). They finally tend to the asymptotic values 1 and as obtained 
in Eqs. (62). 

Let us note that the behavior of the two models, although similar in the limiting case of very high temperatures, 
display some qualitative differences. As can be seen by comparing Eqs. (45) and Eqs. (60), in the Ai model the 
density- pressure ratio does not rise above the usual value 3, contrary to what happens in A 2 . 

We have here a first manifestation of the domination of bosons in this context: their density-pressure ratio goes 
like a constant while the one corresponding to fermions vanishes. Their entropy is also the only one to be considered 
at scales much higher than T cr as they correspond to n = 2 in Eq.(47) while for fermions one has n = 3(see Eqs. (61)). 
Although the ratio we considered is much higher for bosons, one still can not infer this to be the case for each of its 
parts. This will be established for radiation in the next section. 



IV. THE BLACK BODY - RADIATION 



Let us now analyze how radiation gets affected. We shall restrict ourselves to the A2 model for the two reasons 
explained above. As the preceding section showed that the equation of state, for example, presents the same features 
at very high temperatures in the two models, we expect the same to occur here. 

The formula we will obtain in this section will be more relevant to the very early cosmology. The reason is that at 
very high temperatures, the spontaneous symmetry breaking which gives masses to particles has not taken place yet. 
One has thus to consider massless particles, i.e particles with zero chemical potential. 

Once again, we follow the notations of [35]. The important quantities are 



qbo 



^log(l-exp(--^)) =-logZ bo , g/e = -^log(l + exp(--^)) =-logZ /e 



(63) 



where Z is the grand partition function. The entropy is given by 

with <S> = kTlogZ 



8<P 
~ ~~dT 



while the energy and the particle number read 



C/ = E 



r ex P (w) - 1 



7 1 exp (w) - 1 



(64) 



(65) 



It should be noted at this point that considering the A\ model for example, if one did not impose an ultraviolet cut 
off, the periodic dependence of the energy (Eq. (41)) would have led to a divergent energy density. 

For bosons in the A 2 model, the quantity q linked to the partition function by Eq.(63) can be written as 



AttV 



/kT\ 3 f^ 2 * 
V W Jo 



dxx log 



1 — exp — 



1 + 2 TZ x2 , 



(66) 
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while the energy assumes the following form 

n4 /-V5?F 



U = AttV 



(kry 

(he) 3 



f 

Jo 



dx ■ 



1 + 



2 



exp 



1 I 2 ^2 *^ 



(67) 



The particle number admits a similar integral expression. 

For temperatures greater than or comparable to T cr , a Taylor expansion to fourth order leads to the following 
expressions: 



Pbo = &T C 
Pbo = cT c 



2 r-T c r-T /8 , T 112 4 , „ 



8 _T 4 ^T c 

-2+ -V2 1 V2— 

3 T r 15 T 



- -y/2 + 6 — 

3 T 3 T r 



Sb = (J 



where a 



^(^) +^V2^601og- + 116-301og2 



/i 3 c 3 



(68) 



The corresponding quantities for fermions (except the pressure) are dominated by constants: 



-V^log2--V^f^ 

3 s 5 It 



Sf e = a 

The behavior of the energy is depicted in Fig2. 



12 - 
10 - 

8 - 

6 

4 - 

2 



(69) 




0.5 1 1.5 2 2.5 3 
FIG. 2. The energy densities for fermions and bosons are plotted in terms of the temperature. The units adopted are aT cr 
and tS—. 

1 cr 

At temperatures below T cr , the energy density is polynomial (<~ T 4 ). Above T cr , it becomes linear as obtained in 
Eq.(68) for bosons while it goes to a constant for fermions, as shown in Eq.(69). 

The difference between bosons and fermions in the unmodified theory is encoded in the factor 7/8, for the energy 
contributions for example. One sees this is dramatically changed here. In the usual theory, the ultra relativistic gas 
shares the same equation of state with the black body radiation; this feature is also lost when a minimal uncertainty 
in length sets in. 



V. CONCLUSIONS 



We have studied the thermodynamics induced by a non local theory which exhibits a minimal uncertainty in length. 
We have obtained that a new behavior sets in at very high temperatures. The difference between fermions and bosons 
is more important than in the usual case. 

It is worth mentioning some aspects which have not been raised in this work. At the fundamental level, one can 
ask if the concept of spin is relevant in these theories and, in the case the answer is positive, one still has to study 
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the relation between spin and statistics in the new context. As the spin of a particle is defined, in the modern 
approach, through the behavior of its wave function under the Lorentz group, one has to find its generalization in the 
new context. For example, taking the ultimate structure of space-time to be given by a particular noncommutative 
geometry, the relevant algebra is not the Poincare algebra but its q deformation. A notion of spin has been defined in 
these theories and the wave equations for particles of spin 0, 1/2 and 1 have been found [36]. Although the question 
has not been addressed in K.M.M. theory, we hope a similar situation to occur. We expect a generalization of the 
notion of spin which conserves the spin-statistic theorem. 
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